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Abstract. Using deep VLT/ISAAC near- infrared imaging data, we derive the iG-band luminosity function (LF) 
of the z = 1.237 massive X-ray luminous cluster of galaxies RDCS J1252.9— 2927. Photometric redshifts, derived 
from deep multi-wavelength BVRIzJK s data, and calibrated using a large subset of galaxies with spectroscopic 
redshifts, are used to separate the cluster galaxy population from the foreground and background field galaxy 
population. This allows for a simultaneous determination of the characteristic magnitude K* and faint end slope 
a of the LF without having to make an uncertain statistical background subtraction. The derived LF is well 
represented by the Schechter function with K* — 18.54^'^ and a = — 0.64j^'2g. The shape of the bright end of 
the derived LF is similar to that measured at similar restframe wavelengths (in the z-band) in local clusters, but 
the characteristic magnitude is brighter by AM* = —1.40^_q' BS magnitudes. The derived a is similar to the value 
derived in the _ff s -band in the z — 1 cluster of galaxies MG2016+112 but is shallower (at the 2.2<r level) than the 
value measured at similar restframe wavelength (in the z-band) in clusters in the local universe. The brightening 
of the characteristic magnitude and lack of evolution in the shape of the bright end of the LF suggests that the 
massive cluster ellipticals that dominate the bright end of the LF were already in place at z = 1.237, while the 
flattening of the faint end slope suggest that clusters at * ~ 1 contains relatively smaller fractions of low mass 
galaxies than clusters in the local universe. The results presented in this paper are a challenge for semi analytical 
hierarchical models of galaxy formation which predict the characteristic magnitude to grow fainter and the faint 
end slope to steepen with redshift as the massive galaxies break up into their progenitors. The observed evolution 
is consistent with a scenario in which clusters are composed of a population of massive galaxies which formed 
at high redshift [z 3> 1) and subsequently evolved passively, and a population of lower mass galaxies which are 
gradually accreted from the field, primarily at lower redshift (z < 1). 
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1. Introduction 

The galaxy luminosity function (LF) is a powerful tool 
for constraining models of galaxy formation and evolu- 
tion. The LF of both field and cluster galaxies have been 
shown to be well represented by the Schechter function 
(f>(m)dm = n* [l0 ' 4 ( m *- m )] Q+1 e- 10 ° 4(m *" m, dm, where 
to* is the characteristic magnitude of the distribution, a 
is the faint end slope, and 77* is a normalization constant 
describing the space density of galaxies. 

The variation of the LF parameters with galaxy type 
(early- type/late- type), wavelength (UV, optical, near- 
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infrared [NIR] ) and environment (cluster versus field) and 
their evolution with redshift depends on the details of 
mass assembly and star formation and therefore provide 
strong constraints on galaxy evolution models. 

At optical wavelengths the LF of field galaxies and 
cluster galaxies have been derived from a number of stud- 
iGarilli et al.|l999T:|B-lanton et al.l200lUGoto et all 
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de Prooris et all l2003|) . The resulting LF parameters 



varies somewhat from study to study, most likely due to 
differences in the data quality, the employed methods for 
selecting complete samples galaxies and deriving their lu- 
minosity functions. 

The most exten sive studies are d erived from the Sloan 
Digital Sky Survey ijYork et al . 2000). Initial results based 
on the SPSS co mmissioning data ijBlanton et al.l EoOlt 
iGoto et all 12002) disagreed with results from the litera- 
ture, probably due to problems with the photometry, k- 
corrections and evolutionary corrections. Studies based on 
more recent data, however, agree better wit h the literature 
and do not seem to have such problems ijBlanton et, all 
l2003t IPopesso et allEool . These studies show that the 
LF of field galaxies and cluster galaxies are very similar. 
The LF in both cases has brighter characteristic magni- 
tudes and steeper slopes in the redder bands than in the 
bluer bands. The faint end slope is slightly steeper for 
cluster galaxies than for field galaxies in all bands, while 
the characteristic magnitude from these studies appear 
to depend very little on environment. Rich evolved clus- 
ters with cD galaxies tend to have flatter faint end slopes 
than poorer clusters, indicating a deficiency of star form- 
ing dwarf galaxies in rich c l usters llLooez-Cruz et" alll997l 
iDriver fc De Proprisl20oH) . iBarkhouse et all (|2002^ find a 
tendency for the faint end slope, of the LF of local clus- 
ters, to steepen with cluster centric distance, consistent 
with the hypothesis that a large fraction of dwarf galaxies 
near cluster centers are being tidally disrupted. 

While the variation in the LF parameters at optical 
wavelengths are sensitive to the star formation properties 
and morphological mix of the underlying galaxy popula- 
tion, the NIR LF is an excellent probe of the mass function 
of clusters. The if-ban d flux of a galaxy is a good tracer 
of its stellar mass (e.g. iGavazzi et a,l.lfl996h it is not very 
sensitive to star formation and attenuation by dust (com- 
pared to the optical band) and the k-corrections are fairly 
small and relativel y independent of gala xy type even at 
high redshifts (e.g. iMannucci et al.ll200lh . The evolution 
of the mass function is directly predicted from galaxy for- 
mation and evolution models, making the if-band LF of 
field and cluster galaxies very useful for fundamental test- 
ing of the models. 

Locally, the field galaxy if-band LF has been de- 
rived using the 2 Micron All Sky Survey (2MASS, 
IJarrett et all2000fr . Early-type and late- type field galaxies 
have similarly shaped LF but the early-types are brighter 
(more massive) and l ess numerous than the late-types 
( Kochanek et al.ll200lh . The evolution of the field galaxy 
if-band LF has been estimated in three broad redshift 



bins centered on z = 0.5 , 1.0 and 1.5 using the K20 sur- 
vey l|Pozzetti et al.1 l20"03j) . No evolution is found in a to 
z = 1.0, in the highest redshift bin the data are not deep 
enough to constrain it. The characteristic if-band mag- 
nitude of the galaxies in the z — 1 bin (which includes 
galaxies in the redshift range 0.75 < z < 1.3) is consistent 
with a mild luminosity evolution in the if-band to z ~ 1 
of AMk = -0.54 ± 0.12. At z = 1, the bright end of the 
LF is still dominated by ellipticals, and only a small de- 
crease in their space density is observed, indicating that 
the field ellipticals are largely in place at z = 1. 

In a deep s t udy o f the LF of the Coma cluster, 
Ide Propris et all l)l998|) find a characteristic magnitude 
M* K = -24.02 and a faint end slope of a = -0.78. If they 
restrict the fit to "bright" ga l axies (M < M* + 3) they find 
a = -0.9. Ide Propris et all l|l999l ) studied 38 clusters in 
the redshift range 0.1 < z < 0.9 to look for an evolution in 
the if -band cluster LF. Their data was not deep enough to 
constrain the faint end slope , so they fixed it on the Coma 
cluster value a = —0.9, and concentrated on searching for 
an evolution in the properties of the bright end. Using 
the same method as Ide Propris et al ] lll999h the sample 
has subsequently been extended to inc lude a few clusters 
in the redshift range z = 1.0 — 1.2 ijKodama fc Bower] 
l2003tlNakata et all200ll . The conclusion from these stud- 
ies is that the shape of the bright end of the if-band LF 
appears unchanged to z — 1.2 and the evolution of the 
characteristic if-band magnitude with redshift K*(z) is 
well described by passively evolving galaxies assembled at 

A high formation redshift is also inferred from 
the evolution of the colour magnit ude (CM) relation 
of cluster early- types to z ~ 1.3 l|Bower et all 1992; 
Aragon-Salamanca et al.lll993llStanfqrd et alJll998l Il997t 
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2004J). The simplest explanation for the evolution of 



the CM relation is the monolithic collapse models (e.g. 
lEggen et al.lll962^1 in which the bulk of stars are formed 
in a single short duration burst at high redshift and sub- 
sequently evolve passively, but it can be equally well 
explained by semi analyti cal hierarchical models (e.g. 
iKaiiffrnann fc Charlotll99^ in which the stars are formed 
in the disks of of late-type progenitors which later merge 
to form the elliptical galaxies, provided that the merging 
does not trigger significant star formation. 

The if -band LF is predicted to evolve quite differently 
in the two formation scenarios, due to its close relationship 
with the evolution of the mass function. 

If the galaxies were assembled at high redshift and 
subsequently evolved passively, as predicted by the mono- 
lithic models, we expect the characteristic magnitude to 
brighten with redshift (as a consequence of the passive 
evolution of the stars) but the shape of the LF to remain 
the same. If the galaxies were build up through continuous 
merging over a broader redshift interval as predicted in 
hierarchical models, we expect the shape of the LF to 
change with redshift to reflect the accretion history of the 
cluster and the merging history of the cluster galaxies. A 
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deficiency of the brightest galaxies should appear as the 
the most massive galaxies break up into their progenitors, 
resulting in a steepening of the bright end of the LF with 
redshift. If small structures form first and subsequently 
merge to form larger structures, the faint end slope should 
steepen with redshift as galaxies break up into their pro- 
genitors. In many hierarchical models however, the cluster 
population grow by accretion of surrounding structures. 
The accretion history, and the mass spectrum of the 
accreted structures could influ ence the evolution of the 
LF. The Butcher-Oemler effect l|Butcher fc Oemlerlll978|) 
in intermediate redshift clusters is probably due to star 
formation in low mass galaxies which have been accreted 
from the fiel d, induced by intera c tion w it h the cluster en- 
viron m ent llButcher fc Oemlerl Il984t, I Fabrica nt et alJ 
19911 ISmail et all 1998t Ivan Dokkum et all l200oF 
de Propris et al.l \200c ). Such an accretion of primar- 



ily low mass galaxies could cause the faint end slope to 
decline with redshift, as a smaller number of low mass 
galaxies would have had time to be accreted onto the 
cluster environment at high redshift. 

The apparent passive evolution of the bright end of 
the cluster LF to z ~ 1 is a challenge for hierarchical 
models which predict that the characteristic mass should 
be a factor of three sm aller at z ~ 1 than at present 
i Kodama fc Boweill2(~)03f) . None of the studies mentioned 
above have been able to constrain the evolution of a. 
The first constra ints on a at high redshift was derived by 
iToft et alJ ( 2003) who simultaneously fitted K* and a in 
the z = 1 cluster MG2016+112. The derived constraints 
on a arc not strong (a = — 0.601q 33), however it is im- 
portant to keep it free in the fit since it is coupled to the 
derived value of K* and its uncertainties. Furthermore it 
is noted that since a depends on wavelength, and the K- 
band corresponds roughly to restframe z-band at z ~ 1, 
the derived value should be compared to the local value 
measured in the z-band rather than in the iVband when 
studying evolutionary patterns. 

In this paper, we derive firm constraints on K* and a 
for the iVband LF of the massive, X-ray luminous cl us- 
ter of galaxies RDCS J1252.9-2927 fctosati et al|2004) at 
z = 1.237. We apply the method of IToft et all <)2003[) to 
high quality multi-wavelength NIR and optical photomet- 
ric data and extensive spectroscopic follow up observations 
covering a 4' x 4' field around the cluster. 

The outline of this paper is as follows: In [0 we give 
a brief introduction to the data and describe how cluster 
galaxies are separated from the field galaxy population 
using photometric redshifts, calibrated using the spectro- 
scopic data. In fJ3|we derive the LF of the cluster galaxies, 
in 21 we compare to results from studies of lower redshift 
clusters to look for an evolution and compare it to model 
predictions, and in S[5] we summarize and discuss the re- 
sults. Throughout the paper we assume a flat cosmology, 
with Q m = 0.3, Q A = 0.7 and H = 70 kms" 1 Mpc -1 . 
Magnitudes are in the the Vega system. 




K [mag] 

Fig. 1. Binned representation (bin size 0.2 mag) of the 
number of detected objects as a function of i^ s -band 
magnitude. The data are estimated to be complete to 
K, = 22.5. 



2. Building a complete sample of cluster galaxies 

The observational basis of this study consist of high qual- 
ity multi-wavelength NIR and optical photometric data of 
a 4' x 4' field around RDCS J1252.9-2927, supplemented 
by spectroscopic follow up observations, all obtained at 
the ESO Very Large Telescope (with the ISAAC, FORS1 
and FORS2 instruments). The central part of the field 
was imaged for 24 hours in each of the NIR wavebands: 
J s and K s , and for approximately 1 hour in each of the 
optical wavebands: B, V, R, I and z. The full data set 
is described in detail in t hree companion pap ers: a paper 
describing the NIR data (iLidmanet^lJ 2004) a paper de- 
scribing the optical data ijRosati e t al.||20 04[) and a pap er 
describing the spectroscopic data ijDemarco et al.ll2004)) . 

To derive the luminosity function of the cluster galax- 
ies in RDCS J1252. 9-2927 we need a catalog of cluster 
galaxies which is complete to a certain magnitude. In 
Fig. ^ we plot a binned representation of the number of 
detec ted objects in the i f. q -band as function of magni- 
tude l|Lidman et alJEo04l) . From visual inspection of the 
turnover magnitude in this figure we estimate that the K s - 
band data are complete to K s = 22.5. To build a sample 
of cluster galaxies which is complete to this magnitude, 
the redshift distribution of all galaxies in the field must 
be determined in order to separate the cluster galaxy pop- 
ulation from the foreground and background populations 
of galaxies. 

At faint magnitudes it is not feasible to obtain this 
solely through spectroscopic observations. Instead we take 
advantage of the extensive optical and NIR imaging data 
to derive photometric redshifts and calibrate and test 
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these against a large subset of galaxies with spectroscopic 
redshifts. 

2.1. Photometry and photometric redshifts 

We used SExtr actor l|Bertin fc Arnoutslll996[) for object 
detection in each of the available wave bands, for comput- 
ing magnitudes in apertures of 2" and for cross-correlation 
of the resulting catalogues. The (3ct) limiting (Vega) mag- 
nitudes achieved for the imaging in 2" apertures are: 
B = 26.7, V = 26.8, R = 26.0, / = 25.7, z = 24.3, 
J s = 23.9, K s = 22.5. 

To obtain reliable colours we must take into account 
that the imaging in the different wavebands was obtained 
under different seeing conditions. To derive and apply 
these corrections we did the following: first, we com- 
puted magnitudes in apertures of 2" in each of the avail- 
able ("original") waveband images. We then degraded the 
point spread function (PSF) of the original images to 
match the worst seeing condition and recomputed the 2" 
aperture magnitudes in the "degraded" images. The cor- 
rections to be applied to the "original" magnitudes was 
then derived by comparing the magnitudes derived for 
bright stellar objects in the "original" and the "degraded" 
images. 

We used the B avesian photometric redshift (bpz) code 
of lBemte3 ()2000|) to derive photometric redshifts. The ad- 
vantage of the Bayesian approach is the use of a priori 
probabilities by which it is possible to include relevant 
knowledge, such as the expected shape of redshift distri- 
butions and the galaxy type fractions, which can be read- 
ily obtained from existing surveys. For our study we used 
the sa me se t of templates a s tho se described in iBemtezl 
l|2000|) : four lColeman et alJ l|l980|) templates (E/SO, Sbc, 
Scd and Irr) and th e spectra of two starburst galaxies in 
iKinnev et alJ l|l996h . We derived two interpolated SEDs 
between each pair of thes e spec tral types. We used the 
priors derived by iBemted l)2000() from the Hubble deep 
field north (HDFN) and Canada France redshift survey 
(CFRS) catalogues. 

To quantify the reliability of the photometric redshift 
estimation we used a sample of secure sp ectroscopic red- 
shifts from our campaign in this area ijDemarco et alJ 
In the region with multi-wavelength coverage de- 
scribed above, there are 120 such spectroscopically iden- 
tified sources. 

In Fig. we plot z p h ot versus z spec for the spectro- 
scopic sample. The photometric redshifts in general repro- 
duce the spectroscopic redshift well but are systematically 
smaller. This is illustrated in Fig. where we plot the dis- 
tribution of deviations of the photometric redshifts from 
the spectroscopic redshifts. The (sigma clipped) mean of 
the distribution is {z p hot — z spec ) = —0.04 and the stan- 
dard deviation is o = 0.18. If we restrict the analysis to 
the spectroscopic cluster members (the filled histogram 
in Fig. |2J the mean deviation is larger (z p hot — z spec ) = 
—0.13 and the standard deviation is a = 0.16. In Fig§] 



2.0 1 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 r 



0.5 - 



o.o i l_j i i i i i i i i i i i i i 

0.0 0.5 1.0 1.5 2.0 

zspec 

Fig. 2. Photometric redshift (z p hot) versus spectro- 
scopic redshift (z spec ) for the full spectroscopic sample. 
Immediately recognizable is the cluster of galaxies at 
z = 1.237. 
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Fig. 3. Distribution of deviations of the photometric red- 
shifts from the spectroscopic redshifts for the spectro- 
scopic cluster member galaxies (filled histogram) and for 
the full spectroscopic sample (open histogram). 



we plot (zphot — z spec) as a function of redshift to investi- 
gate whether this is a consequence of a redshift dependent 
systematic photometric redshift error, which could be in- 
troduced by small errors in the photometric zero-points of 
one or more of the filters. No significant dependency on 
redshift is observed. In the following we add 0.04 to all the 
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Fig. 4. Mean deviation of the derived photometric redshift 
from the spectroscopic redshift as a function of spectro- 
scopic redshift, calculated in S — 0.1 bins. No significant 
dependency on redshift is observed. The dotted line marks 
the mean offset —0.04 of the full sample. 

derived redshift to empirically correct for the systematic 
deviation. 

To compare the photometric redshifts errors with those 
achieved in the literature we define for each object the 
reduced error in the photometric redshift estimation as 
S x = {zphot - z spec )/(l + z spec ). For the full sample we 
achieve a mean offset (S z }—0.02 and an standard devi- 
ation of er((5 2 )=0.11. If we restrict our analysis to the 
cluster members we have (S Z )—0.Q5 and a(5 z )=0.09. The 
accuracy of our photomet ric redshifts are com parable to 
the accuracy achieved by iBarger et all 11200.1) m a simi- 
lar Bayesian photometric redshift study of galaxies in the 
Chandra deep field north. For the HDFN, iBem'te d l|2000|) 
achieve a higher accuracy: a(5 z ) — 0.06, probably due to 
higher precession Hubble Space Telescope (HST) photom- 
etry. 

2.2. The photometric cluster member sample 

In this section we describe how the cluster galaxy popula- 
tion is separated from the foreground and background field 
galaxy population using the photometric redshift. We take 
into account the uncertainty in the photometric redshift 
estimation by allowing galaxies with photometric redshifts 
within Az of the cluster redshift z c i to be classified as 
"photometric members". The choice of Az is a trade off 
between cluster galaxy completeness and field galaxy pol- 
lution. The interval must be sufficiently broad to include 
as many "real" cluster galaxies as possible and sufficiently 
narrow to minimize pollution from field galaxies. From 
Fig. 12141 it can be seen that most cluster member galaxies 



Fig. 5. Completeness function of the photometric mem- 
ber sample, calculated as the fraction of simulated z = 
1.237 cluster galaxies recovered by the photometric red- 
shift analysis as a function of magnitude. 

at z = 1.237 are expected to have photometric redshifts 
in the range \z p hot — z c i\ < 0.3. Since we are aiming at 
building a sample of cluster galaxies which is complete to 
K s = 22.5 we prioritize not to exclude any cluster mem- 
bers and choose Az = 0.3. In the following, galaxies with 
0.935 < z p h < 1.535 are thus classified as photometric 
members. 

To estimate the "completeness function" of the photo- 
metric member sample (the fraction of z — 1.237 cluster 
galaxies recovered by the photometric redshift analysis in 
the Az interval, as a function of magnitude) we did the fol- 
lowing: for each 5m = 0.5 mag bin in the observed range of 
magnitudes K s = 16 — 23, we generated a catalog of 1000 
galaxies at z — 1.237 as they would appear in a dataset 
with the same band-passes and limiting magnitudes as 
the data of RDCS J1252. 9-2927. The catalog was gener- 
ated using the make-ca talog code in the hyperz package 
l|Bolzonella et alJ l2000) and were drawn randomly from 7 
template spectral types (E to Im). We then applied BPZ 
to the catalogs. In Fig. OD we plot the fraction of input 
(z = 1.237) galaxies with derived photometric redshifts 
within 0.935 < z p h < 1.535 as a function of if s -band 
magnitude. Down to a magnitude of K s = 21.5 about 90% 
of the galaxies are recovered. At K s = 22.5 the recovery 
rate is 70%. The reason why BPZ appears to be performing 
slightly worse at the brightest magnitudes, compared to at 
fainter magnitudes, is that the simulated catalogs contain 
sizable fractions of bright late- type galaxies which are very 
rare at z = 1.237, and therefore are assigned small proba- 
bilities by BPZ, resulting in a some cases, in the most likely 
spectrum being a lower redshift spectrum of earlier type. 
The redshift evolution of the field galaxy luminosity func- 
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Fig. 6. Contours of the smoothed density distribution of 
photometric member galaxies overlaid on the ~ 4 ' x 4' 
mosaic i^ s -band ISAAC image (jLidman et al.ll200^ . The 
contours are 2 times (dashed line) and 4-10 times (full line) 
the density of galaxie s in the HDFN with K s < 2 2.5 and 
0.935 < z ph < 1.535 (|Fernandez-Soto et alJll99fih . Inside 
the dotted circle with radius 65" the contamination from 
field galaxies in the photometric member sample is 25%. 
This circle marks the aperture inside which we derive the 
cluster galaxy luminosity function of RDCS J1252.9— 2927 
in S£0 

tion, and its cosmic variance is not known with sufficient 
accuracy to correctly incorporate the effects of pollution 
of field galaxies in the completeness function analysis. We 
return to the issue of field contamination in M2.3I and 

2.3. Spatial distribution of photometric member 
galaxies 

To investigate the spatial distribution of the photometric 
members galaxies we constructed an image with the pixel 
value at the centroid of the photometric cluster members 
equal to one and the remaining pixel values equal to zero. 
This image was then smoothed with a Gaussian kernel 
with FWHM=60 pixels (~ 8"8). In Fig. El we overlay con- 
tours of the smoothed density distribution of the photo- 
metric member galaxies on the K s -baiid image of the clus- 
ter. The contours are 2 — 10 times the density of galaxies 
in the HDFN with K» < 22.5 and 0.935 < z ph < 1.535 



1 Fernandez-Soto et alJll99^. 



There is a large concentration of photometric member 
galaxies in the central part of the field. In this region the 
photometric member sample must thus be dominated by 
cluster galaxies. The density of photometric members in 
the outer parts of the field is comparable to the density 
in the HDFN suggesting that field galaxy pollution could 




Right Ascension 

Fig. 7. Map of the smoothed if s -band light of the photo- 
metric member galaxies. The map was generated by cre- 
ating a version of the K-band image in which pixels inside 
squares (of size 20 pixx20 pix) centered on the photo- 
metric members were left at their original values and the 
remaining pixels were set equal to zero. This image was 
then smoothed with a Gaussian kernel with FWHM=60 
pixels. The circle marks the region inside which the LF is 
derived. 



make a significant contribution to the photometric mem- 
ber sample in these regions. We discuss this further in ^ 
In addition to the number-density distribution of the 
photometric member galaxies we also investigated the dis- 
tribution of their if s -band light. In Fig. [7| we show a 
smoothed map of the if s -band light of the photometric 
member galaxies. Since the i^ s -band flux traces the stellar 
mass of the cluster galaxies, this map can be be thought 
of as a stellar mass map of the cluster, which can be com- 
pared to maps of the gas mass derived from X-ray obser- 
vations, and maps of the dark matter derived from weak 
lensing observations. With this in mind it is interesting to 
note that, in addition to a strong concentration of light 
in the center of the cluster, the light is extended in the 
East- West direction. This extension is a lso seen seen in 
both X-ray and weak-lensing mass maps l)Lombardi et all 
l2004h . 

3. Luminosity function 

We now have a sample of cluster galaxies which is com- 
plete to K s = 22.5, or rather a sample which we know 
how to correct for incompleteness, and we can derive the 
K s -band luminosity function of the cluster galaxies with- 
out having to make uncertain statistical corrections to ac- 
count for foreground and background field galaxy contam- 
ination. The choice of Az in W2. 21 ensures that all spectro- 
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scopically confirmed cluster members are included in the 
photometric member sample. To take full advantage of 
the data we apply the m aximum likelihood technique of 
ISchechter fc Presa l|l976fl directly to the luminosity distri- 
bution of the cluster galaxies rather than doing a "least 
squares" fit to a binned representation. In this way we in- 
clude the information that in many magnitude intervals no 
galaxies are found, and do not make the assumption of the 
X 2 method that the underlying distribution is Gaussian. 
The "incompleteness" of the photometric redshift selec- 
tion is taken into account through the completeness func- 
tion shown in Fig. [S] This in turn leads to more realistic 
error bars. For more details o f the method we refer to the 
appendix of lToft et all l|2003|) . 

In t]2.3l we argued that pollution from field galaxies in 
the photometric member sample could be significant in 
the outer parts of the 4' x 4' field. The amount of pollu- 
tion can be estimated from the subsample of photomet- 
ric member galaxies with spectroscopic redshifts. These 
galaxies are likely to constitute a relatively fair sample of 
galaxies in the relevant redshift range, since the selection 
of the spectroscopic sample was designed not to introduce 
biases on the cluster galaxy populations, while minimiz- 
ing the pollution of field galaxies. This was accomplished 
by targeting galaxies with K s < 21, J — K s < 2.1 and 
R — K s > 3. Such criteria do not penalize cluster galaxies, 
since at z = 1.237 even the latest types are redder than 
R — K s = 3 and early types are bluer than J — K s = 2.1, 
however field contamination is significantly reduced. 

In FiglHl the bottom panel shows the contamination 
derived in apertures with increasing radius (centered on 
the cluster core) . In the central parts of the field the con- 
tamination is modest. The pollution within 45" is ~ 10%, 
within 65" it is - 25%, and within 120" it is 50%. To 
investigate the effects of the pollution on our analysis we 
derived the luminosity function of the photometric mem- 
ber sample in apertures of increasing size. The top and 
middle panel in FiglHl shows the variation of the luminos- 
ity function parameters with aperture radius. 

There is a tendency for K* to be slightly brighter and 
a to be slightly smaller (less negative) in apertures en- 
compassing only the central regions where the contam- 
ination is small, compared to in larger apertures where 
the contamination is more pronounced, but the effect is 
barely significant since the error bars are larger in the 
smaller apertures due to the smaller number of galaxies. 
Part of the effect could be caused by intrinsic variation in 
the properties of the cluster galaxy LF with cluster cen- 
tric distance, and part of it could be a consequence of 
field galaxy pollution. Since the effect is not statistically 
significant however, field galaxy pollution is not likely to 
significantly affect our results. 

Based on Fig|H|we limit our cluster galaxy LF analy- 
sis to galaxies within 65" of the cluster center, in order 
to maximize the number of galaxies while minimizing the 
pollution. The photometric member sample contains 100 
galaxies within this distance, including 19 spectroscopi- 
cally confirmed cluster members and 7 interlopers with 
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Fig. 8. The top and middle panel shows the K* and a 
parameters of the LF of the photometric member galax- 
ies derived in apertures of increasing size, centered on 
the cluster core. Note that the two parameters are not 
independent, so the error bars in the two plots are cor- 
related. The bottom panel shows the contamination in 
apertures of increasing size estimated from the subsample 
of photometric member galaxies with spectroscopic red- 
shifts. The contamination is calculated as the number of 
interlopers, i.e. photometric member galaxies with spec- 
troscopic redshifts different from the cluster redshift (de- 



fined as 1.22 < 



< 1.25), divided by the total number 



of photometric member galaxies with spectroscopic red- 
shifts. Inside the aperture with radius 65" (marked by the 
vertical dashed line), the contamination is 25%. 



spectroscopic redshifts different from z c \. We remove the 
known interlopers from the sample and study the lumi- 
nosity function of the remaining 93 galaxies in detail in 
the following. 

The Schechter function provides a good fit to the data. 
In Fig. El we plot 1 — 3a likelihood contours of the two 
Schechter function parameters. There is some degeneracy 
between the two parameters, but we are able to put firm 
constraints on both the characteristic magnitude K* = 
18.54±3:|| and the faint end slope a = -0.64tg;|£. Such 
accuracy is unprecedented at these redshifts. 

Since our analysis is based on photometric redshifts 
which can be sensitive to photometric errors, the results 
could potentially be affected by small changes in the pho- 
tometry. To find out whether this is the case, we carried 
out a series of Monte Carlo simulations to investigate how 
sensitive the photometric redshift distribution is to the 
photometric errors. We generated 10 realizations of the full 
dataset, by randomly perturbing the photometry of the 
galaxies within their la error bars. We then derived their 
photometric redshift distribution, defined a photometric 
cluster member sample, removed known interlopers and 
derived their LF in exactly the same way as for the original 
data set. The best fitting LF parameters of the perturbed 
datasets (represented by small symbols in FigEJl all fall 
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within the la contour of the original dataset, indicating 
that the LF analysis is robust with respect to photometric 
perturbations, and that the effect of field galaxies scatter- 
ing in and out of the Az interval as the photometry is per- 
turbed is small, otherwise we would expect to see larger 
variations in the LF parameters. 

In Fig. EH we plot the best fitting Schechter function 
and a binned representation of the data. For comparison, 
we plot the local z-band cluster galaxy LF, rather than 
the local A-band cluster galaxy LF, since the observed 
A s -band corresponds roughly to restframe z-band at 
z = 1. 237. Following the method of lvan Dokkum fc Franxl 
l)l996|) absolute z-band magnitudes M z can be related to 
the observed A s -band magnitudes through: 

M zAB = K s>AB -5*log(d L /W)+2.5log(l+z)+l3(H-K) AB , 

where z = 1.237 is the redshift, dj, is the luminosity dis- 
tance in parsecs and (3{H — K)ab is a colour term to com- 
pensate for the fact that the redshifted z-band does not 
match the observed A s -band exactly. The basic assump- 
tion made to derive this expression is that the flux at the 
redshifted z-band can be related to the observed H and 
A s -band flux by F v {u Zband {z)) = F v {v H f F v {v K f-^ . We 
adopt the value (H — K)ab{z = 1. 25) = 0.57 predicted 
by th e passive evolution models of iKodama fc Arimotol 
l)l997|) for an elliptical L* galaxy formed at 2/ = 3 
(the mean formation re dshift of cluster ellip ticals derived 
from their CM relation. iLidman et alJl2004l) and calculate 
— 0.35 by assuming a simple power-law for the shape of 
the SED between uh and vk- 

Contamination from faint field galaxies is not likely 
to significantly affect the results derived for the cluster 
galaxy LF in the central part of the field, where the pho- 
tometric member sample is dominated by bright galaxies, 
but we note that the main effect of such an contamination 
would be to overestimate the faint end slope, making the 
derived faint end slope a formal upper limit to the intrinsic 
faint end slope of the cluster galaxy LF. 

4. LF evolution 

The observed As-band corresponds roughly to restframe 
z-band at z = 1.237, so the evolution of the z-band LF 
can be constrained by comparing the observed A s -band 
LF with the local z-band cluster LF. This is done in Fig 1 101 
where we plot the derived A s -band LF and the z-band LF 
of clusters in the local universe. For the local LF wc adopt 
the param eters M* = -22.31 ±0 .20 and a = -1.14±0.08 
derived bv lPopesso et alJ l)2004(l for the composite LF of 
X-ray clusters in the Sloan Digital Sky Survey (SDSS). 
From Fig EH it can be seen that the characteristic magni- 
tude of the restframe z-band LF at z = 1.237 is brighter 
than at z = 0. By transforming the observed A* to M* us- 
ing Equation^an evolution of AM Z = — 1.40j^58 ma S is 
derived for the restframe z-band characteristic magnitude. 
The observed A"* can be converted to absolute (restframe) 
A s -band magnitude by applying a k-correction: 

M Ks = K s - 5 * log(d L /10) - k Ka {z = 1.237). (2) 



20 



(1) 




Fig. 9. Contour plot showing the constraints on the 
Schechter function parameters derived from the maximum 
likelihood analysis. The big cross marks the best fitting 
parameters, and the curves represent 1 — 3ct confidence 
levels. The small crosses marks the best fitting parame- 
ters for the catalogs with perturbed photometry (see text). 
The dashed line mark the la upper limit for the faint end 
slope a — —1.14 ± 0.08 derived in lo cal clusters at similar 
restframe wavelength (in the z-band. Ipopesso et alJl2004|) 

If we ado pt the UkAz = 1.25) = —0.68 for the k- 
correction ijMannucci et alJ l200l|) and compare to M£ 
derived locally in the Coma cluster, we derive an evolu- 
tion in the restframe A-band characteristic magnitude of 
AM K = -1.4 ± 0.6 mag (assuming that M* K = M£J, 
in agreement with the evolution derived in the restframe 
z-band. The errors quoted for AM Z and AMk are domi- 
nated by the uncertainty in the derived A*. In Fie llll we 
compare the derived A* with values derived for clusters 
at lower redshift and the predictions of passive evolution 
models f or the evolution of an L* ga laxy formed at Zf = 2, 
3 and 5 ijKodama fc Arimotolll997|) . The observed evolu- 
tion of A* to z = 1.237 is consistent with what is expected 
for a passively evolving population of galaxies formed at 
z f Z2. 

The derived faint end slope a = — 0.64+q'25 is similar 
to the value a = — 0. 60^33 deri ved in the JC-ba nd in 
the z = 1 cluster MG2016+112 (jToft et alJEoO^) . but 
shallower (at the 2.2<r level, see Fig|5J than the value 
a = — 1.14±0.08 derived at simil ar restframe wavelen gths 
(in the z-band) in local clusters IjPonesso et al J I2004) . in- 
dicating that clusters at z ~ 1 contains relatively smaller 
fractions of low mass galaxies than clusters in the local 
universe. This is in disagreement with semi analytical hi- 
erarchical models which predict the faint end slope to 
steepen with redshift as the massive cluster galaxies break 
up into their progenitors. 
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Fig. 10. The full curve is the _ftT s -band LF of the pho- 
tometric member galaxies within 65" of the cluster cen- 
ter, represented by the best fitting Schechter function, 
with parameters K* = 18.54±°;^ and a = -0.64^;^. 
The dashed curve is the local z-band cluster galaxy LF 
which has M* = - 22.31 ± 0.20 and a = -1.14 ± 0.08 
ijPopesso et aljEfJO^) . The filled symbols are a binned rep- 
resentation of the raw counts, while the open symbols have 
been corrected for the incompleteness of the photomet- 
ric redshift selection, using the completeness function in 
Fig El The conversion between observed if s -band magni- 
tude and restframe z-band magnitude is given in the text. 



5. Summary and discussion 

In this paper, we have taken advantage of an extensive 
NIR and optical dataset of the massive, X-ray luminous 
cluster of galaxies RDCS J1252.9-2927 to derive the first 
secure constraints on the shape of the X s -band LF at 
z > 1. The LF was found to be well represented by 
the Schechter function over the observed range of cluster 
galaxy magnitudes: K s — [17.0 — 22.5]. 

We tested our analysis for the influence of photometric 
errors and pollution from field galaxies and found our re- 
sults to be robust and relatively insensitive to the effects 
of field galaxy pollution in the central parts of the field 
where the constraints of the cluster galaxy LF is derived. 



The characteristic magnitude K* 



18.54 



+0.45 
-0.55 



AM Z = lAOtoil magnitudes brighter than the character- 
istic magnitude measured for clusters in the local universe 
at similar restframe wavelengths (in the z-band). 

This is consistent with studies of the fundamental 
plane in a cluster at similar redshift where a luminos- 
ity evolution of A Mr = —1.50 ± 0.13 was found in the 
restframe i?-band II van Dokkum fe Stanforc3ll200l . 

Apart from being shifted to systematically brighter 
magnitudes, the shape of the bright end of the LF at 



Fig. 11. Evolution of K* with redshift. The curves rep- 
resent the evolution of an L* galaxy formed at Zf = 
2, 3 a nd 5 predicted by t he pa ssive evolution mod- 
els of iKodama fc Arimotcl l|l997j) . no rmalized to the 
Com a cluster which have K* = 10.9 ijde Propris et alJ 
Il998j) . The "no evoluti on" predictions are cal culated us- 
ing the k-corrections of iMannucci et ahl l)200l|) . The data 
points from t he literature have been derived with a fixed 
-0.9. jde Propris et alJ Il999t {Nakata et all l200li 



lEllis fc JonedboO^lKodama fc Bowerl2003|) except for the 
z = 1 value which was was derived with a as a free param- 
eter l|Toft et al Jl200^ . The z = 1 data point marked by a 
star, is the mean K* of 5 clusters in the redshift range 
08< z < 1.1 llKoda.ma fc Bower! IS IemTs fc .Toned 
2004), with error bars added in quadrature. 



z = 1.237 appears similar to in the local universe. Since 
the X-band LF is a good tracer of the stellar mass function 
of the cluster galaxies, this suggests that the massive ellip- 
tical that dominate the bright end of the LF were already 
in place at z = 1.237. This is a challenge for hierarchi- 
cal models which predict the bright end of the i^-band 
LF to steepen and K* to become fainter at high redshift 
as the massive galaxies break up into their progenitors. 
At z ~ 1, current hierarchical models predict the charac- 
teristic mass (closely related to the characteristic if-band 
magnitude) to be a t hird of that in the local universe 
l|Kodama fc Bowerll2003|) , which is clearly not the case for 
RDCS J1252.9-2927. 

The brightening of the characteristic magnitude, and 
lack of evolution in the shape of the bright end of the LF 
to redshift z = 1.237 is consistent with a simple forma- 
tion scenario in which the massive elliptical galaxies that 
dominate the bright end of the if-band LF are passively 
evolving systems assembled at high redshift Zf ~ 3. 

This formation scenario is also in agreement with the 
observed properties of the CM-relation of elliptical galax- 
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ies in RDCS J1252.9-2927, which is identical to the CM 
relation found in local clusters in terms of slope and scat- 
ter, but bluer on average, consistent with old populations 
of stars formed at 2 .7 < z f < 3.6 (jLidman et alJ 12004 
Blakcslec et al 1l2003l) . From the evolution of the CM rela- 
tionship alone it is not possible to distinguish between for- 
mation scenarios where the old stars are formed in mono- 
lithic collapse of the elliptical galaxies at high redshift, 
and scenarios where they are formed in the disks of less 
massive late-type galaxies which later merge to form the 
ellipticals, as long as the merging does not trigger signif- 
icant star formation. From the lack of evolution in the 
shape of the bright end of the -ftf-band LF we can however 
deduce that if the massive ellipticals in clusters formed 
through merging, it took place at higher redshifts (z » 1) 
than is predicted by current semi analytical models. 

The results derived here for the evolution of the bright 
end of the cluster galaxy iCj-band LF are similar to the 
results derived from the K20 survey for the evolut ion of 
the field galaxy K s -band LF ((Pozzetti et al.ll2008j) . The 
magnitude of the luminosity evolution (AMk — — 1.4±0.6 
mag to z = 1.237) cannot be directly compared to the 
value derived in the field (AM K = -0.54 ±0.12 mag to 
z = l. lPozzetti et al . 2003) since the latter is a mean value 
derived in a broad redshift interval, centered on a lower 
redshift, but we note that they are broadly consistent. 

It is interesting to note that a population of massive, 
evolved and highly clustered galaxies in the redshift range 
2 z ^ 4 have recently been discovered in deep NIR 
observatio ns ( Fran x et al. 2003; Ivan Dokkum et all 2003; 
iDaddi et all2 003). These are good candidates for progeni- 
tors of the old, massive, passively evolving elliptical galax- 
ies observed in z ~ 1 clusters. 

The depth of our data allow us to trace the LF down to 
4 magnitudes below K*, which is unprecedented at these 
redshifts, and put firm constraints on the faint end slope of 
the LF. We derive a slope a = — 0.64^Q25 which is similar 
to the slope measured in the if s -band in the MG2016+112 
cluster at z = 1 but shallower than the value measured at 
similar restframe wavelengths in clusters in the local uni- 
verse. The observed evolution in the faint end slope can be 
interpreted as an evolution in the stellar mass spectrum of 
the low mass galaxies that dominate the faint end of the 
iGband LF, indicating that high redshift clusters contains 
relatively smaller fractions of low mass galaxies than clus- 
ters in the local universe. This results is in disagreement 
with semi analytical hierarchical assembly models, which 
predict the relative fraction of low mass galaxies to in- 
crease, and the faint end slope of the NIR LF to steepen 
with redshift as massive cluster galaxies break up into 
lower mass progenitors. 

The results presented here are consistent with a sce- 
nario in which clusters are composed of a population of 
massive cluster galaxies which were formed at high red- 
shift (z S> 1) and subsequently evolved passively with lit- 
tle additional star formation and interaction, and a pop- 
ulation of lower mass galaxies which are continuously ac- 
creted from the field, primarily at lower redshift (z < 1). 
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